Detection of Binding Sites on SARS-CoV-2 Spike Protein Receptor-Binding Domain by Molecular Dynamics Simulations in Mixed Solvents.

The novel SARS-CoV-2 uses ACE2 (Angiotensin-Converting Enzyme 2) receptor as an entry point. Insights on S protein receptor-binding domain (RBD) interaction with ACE2 receptor and drug repurposing has accelerated drug discovery for the novel SARS-CoV-2 infection. Finding small molecule binding sites in S protein and ACE2 interface is crucial in search of effective drugs to prevent viral entry. In this study, we employed molecular dynamics simulations in mixed solvents together with virtual screening to identify small molecules that could be potential inhibitors of S protein -ACE2 interaction. Observation of organic probe molecule localization during the simulations revealed multiple sites at the S protein surface related to small molecule, antibody, and ACE2 binding. In addition, a novel conformation of the S protein was discovered that could be stabilized by small molecules to inhibit attachment to ACE2. The most promising binding site on RBD-ACE2 interface was targeted with virtual screening and top-ranked compounds (DB08248, DB02651, DB03714, and DB14826) are suggested for experimental testing. The protocol described here offers an extremely fast method for characterizing key proteins of a novel pathogen and for the identification of compounds that could inhibit or accelerate spreading of the disease.

Screening of Natural Products Targeting SARS-CoV-2-ACE2 Receptor Interface - A MixMD Based HTVS Pipeline.

The COVID-19 pandemic, caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a severe global health crisis now. SARS-CoV-2 utilizes its Spike protein receptor-binding domain (S-protein) to invade human cell through binding to Angiotensin-Converting Enzyme 2 receptor (ACE2). S-protein is the key target for many therapeutics and vaccines. Potential S-protein-ACE2 fusion inhibitor is expected to block the virus entry into the host cell. In many countries, traditional practices, based on natural products (NPs) have been in use to slow down COVID-19 infection. In this study, a protocol was applied that combines mixed solvent molecular dynamics simulations (MixMD) with high-throughput virtual screening (HTVS) to search NPs to block SARS-CoV-2 entry into the human cell. MixMD simulations were employed to discover the most promising stable binding conformations of drug-like probes in the S-protein-ACE2 interface. Detected stable sites were used for HTVs of 612093 NPs to identify molecules that could interfere with the S-protein-ACE2 interaction. In total, 19 NPs were selected with rescoring model. These top-ranked NP-S-protein complexes were subjected to classical MD simulations for 300 ns (3 replicates of 100 ns) to estimate the stability and affinity of binding. Three compounds, ZINC000002128789, ZINC000002159944 and SN00059335, showed better stability in all MD runs, of which ZINC000002128789 was predicted to have the highest binding affinity, suggesting that it could be effective modulator in RBD-ACE2 interface to prevent SARS-CoV-2 infection. Our results support that NPs may provide tools to fight COVID-19.